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CN • investigated theoretically. Using an approach based on the Korteveg-de-Vries equation, it is shown 
^ ] that the critical value of the negative ion density that separates the domains of existence of compressi- 
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on and rarefaction solitons depends continuously on the dust density. A modified Korteveg-de Vries 

■ equation for the critical density is derived in the higher order of the expansion in the small parameter. 

'p_i! It is found that the nonlinear coefficient of this equation is positive for any values of the dust density 

O • and the masses of positive and negative ions. For the case where the negative ion density is close to 

^ ■ its critical value, a soliton solution is found that takes into account both the quadratic and cubic 

. nonlinearities. The propagation of a solitary wave of arbitrary amplitude is investigated by the quasi- 

^ I potential method. It is shown that the range of the dust densities around the critical value within 

^ ' which solitary waves with positive and negative potentials can exist simultaneously is relatively wide. 
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PACS: 52.35. -g, 52.35. -Fp 

1 INTRODUCTION 

It is well known that the presence of negative ions in plasma substantially modifies ion-acoustic 
solitons described by the Korteveg-de-Vries (KdV) equation. This is related to the fact that the 
nonlinear coefficient of the equation can vanish at a certain critical value of the negative ions density. 
If this density is below the critical value, then the solution has the form of a solitary compression 
wave; otherwise, it has the form of a rarefaction wave. To describe the nonlinear properties of an 
ion-acoustic wave at negative ion densities close to the critical density, it is necessary to take into 
account the higher orders of the expansion in the small parameter. In doing so, Watanabe [1] arrived 
at a modified KdV (MKdV) equation, which was then analyzed in [2]-[5]. Mishra and Chhabra [6] 
showed that, in a plasma with a finite ion temperature, there can exist two types of ion-acoustic 
waves— fast and slow ones. In [7], the influence of nonisothermal electrons on the properties of ion- 
acoustic solitons in plasmas containing positive and negative ions was studied by the quasi-potential 

1 



method. The propagation of nonlinear compression and rarefaction ion-acoustic waves in a bi-ion 
plasma was studied by McKenzy et al. [8]-[9] in a gas-dynamic approach, which is alternative to the 
quasi-potential method. 

In the present work, the propagation of ion-acoustic solitons in a dusty plasma containing both 
positive and negative ions is investigated theoretically. It is well known [10] that, when a plasma 
with one positive ion species contains a negatively charged dust, there is a critical dust density below 
(above) which an ion-acoustic soliton has the form of a compression (rarefaction) wave. When the 
plasma also contains negative ions, the problem has two independent parameters— the negative ion 
density and the density of dust, which, depending on the external conditions, can be charged ether 
negatively or positively. As a result, at a given dust density, the critical value of the negative ion 
density can be shifted or can even be absent at all. Moreover, for given masses and temperatures of 
the plasma ion components, the critical value of the negative ion density depends continuously on 
the dust density. Therefore, the nonlinear coefficient in the KdV equation vanishes on a continuous 
line in the „dust density-negative ion density" plane. This line separates the domains of existence of 
compression and rarefaction ion-acoustic solitons. 

The vanishing of the nonlinear coefficient in the KdV equation at the points corresponding to the 
critical density unambiguously relates the dust density to the negative ion density. In a dusty plasma, 
the nonlinear coefficient in the MKdV equation is a function of one of the above parameters (e.g., 
of the dust density) rather than merely a number. Therefore, it is interesting to examine whether 
this coefficient can actually vanish and change its sign. The vanishing of the nonlinear coefficient in 
the MKdV equation would mean the existence of a double critical point at which of both the first- 
and second-order nonlinearities disappear. To find out whether the nonlinear coefficient changes its 
sign is important for understanding the properties of solutions to the MKdV equation. When the 
coefficient is positive, the solution has the form of a soliton. In contrast, when it is negative, a double 
layer can form in the system. The results of our study show that the nonlinear coefficient is positive 
at any values of the dust density and the masses of positive and negative ions. Therefore, in a bi-ion 
dusty plasma, the solution near the critical points has the form of a soliton. 

In studying the propagation of an ion-acoustic soliton in a dusty plasma, we will assume that the 
charge of a dust grain in the soliton field remains constant. This is true if the charging frequency 
Vch — a^^pi/V^^Di (a - is the dust grain radius, ujpi — (47re^nio/?7T.j)^/^ - is the plasma ion frequency, 
and Xdi — {Ti/4:ne'^nioy^^ - is the ion Debye radius) is considerably smaller than the reciprocal of the 
time during which the soliton passes through a given point in the plasma. This time can be estimated 
as Ti ~ L/cs ~ NXDi/cs = Nu^Ti/Te < Nuj^^, where L - is the soliton length and N = L/X^i. 
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The condition for the grain charge to be constant is then satisfied when Nvch/(jJpi ~ Na/Xm <S 1- 
For a ~ 10~^ cm, Tj ~ 1 eV, and N — 100, the above assumption is valid if lO'^nH^ <S 1, which 
corresponds to a fairly low plasma density (rzjo < 10^). In a dense plasma, the charge of dust grains 
in the field of an ion-acoustic soliton cannot be considered constant. In this case, either ion-acoustic 
shocks [11] or weakly-dissipative ion-acoustic dust solitons [12] can form, depending on the electron 
distribution function. 

The present paper is organized as follows. In Section 2, we write out a set of basic dimensionless 
equations. In Section 3, the KdV equation is derived using the method of perturbations and the 
dependence of the amplitude of a solitary wave on the plasma parameters is investigated numerically. 
In Section 4, the MKdV equation is derived, its nonlinear coefficient is analyzed, and a solution is 
obtained that describes an ion-acoustic wave in a plasma at densities close the critical density. In 
Section 5, the propagation of a wave of arbitrary amplitude is investigated by the quasi-potential 
method. Finally, in Section 6, the main results of this study are summarized. 



2 BASIC EQUATIONS 

Let us consider a dusty plasma containing positive and negative ions. For definiteness, we assume 
that the ions are singly charged. The problem to be solved is described by the set of equations 
that includes the equations of continuity and motion for both ion components, as well as Poisson's 
equation: 

^ ^ y^^^ ^ - ^O^m— (2) 
dt dx J dx dx ^ 

_ + _ ^ 0, (3) 

d d\ I (d(j) dn2\ , . 

— = e"^ - i^iui + H2n2 + A*. (5) 

Here, the coordinate x, the time t, and the velocities vi, V2 are normalized to the electron Debye radius 
Xd — (47rneoe^/Te)^/^, the inverse plasma frequency uj~^ = (mi/47rneoe^)^/^, and the ion-acoustic 
velocity = (Te/mi)^/^ respectively; fj,i = nao/^T-eo, = ni2o/neo, (p = e^'/Te, 9i = Tn/T^, 
O2 = Ti2/Te, Q — m2/mi, Te,Tii,Ti2 are the temperatures of electrons and positive and negative 
ions, respectively. The densities of electrons and positive and negative ions are normalized to their 
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equilibrium values, rigo, ^T-iiO) ^120, respectively. We assume that the wave propagates adiabatically, so 
the electrons density is Ug — exp(0). The parameters /ii and ^2 are related by the quasineutrality 
condition fii = 1 + fi2 + 1^-, where /i = Znao/ngo, Z is the charge number of a dust grain, and nao is 
the dust number density. 

3 KORTEVEG-DE-VRIES EQUATION 

We will search for the solution to the above set of equations in the form 

A = Ao + £Ai+£2A2 + £^A3 + ..., (6) 

where A = {ni,vi,n2,V2,(f)), Aq = (1,0,1,0,0). To derive the Korteveg-de-Vries equation, we make 
use of the standard method of coordinate stretching, 

e = r^e^/H. (7) 

In the first order in e, the set of equations is written as 

-F%i + § + 3^,%i = 0, (9) 

-I'%i4§ + 3§%1=0. (11) 

01 - yWinii + fi2n2i = 0. (12) 

Integrating Eqs. (8)-(ll) with the initial conditions Vii{^ — > ±00) = 0, i'2i(^ — > ±00) = 0, nii(^ — > 
±00) = 0, 7121 (^ ±00) = and substituting the expressions for riu and 7121 into Eq. (12), we obtain 
the following expressions for the first-order quantities, as well as the relationship that determines the 
linear velocity of the wave: 

01 V ^ 

01 V . .... 

/^l ^ ^ ^ 
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In the second order in £, we have 



drill ..9ni2 , dvi2 d 



dr 



Or 



dn2i 
dr 



V 



dn 



22 



dv 



22 



d 



(^21^21, 



0, 



dV2l ^.dV22 , 

- + V2l 



dr 



0^ 

de 



dv2i 1 d(j)2 302 dn22 3^2 dn2i 

= 02 + 0l/2 - flini2 + fi2n22- 



0, 



Excluding the term dvi2/d^ from Eqs. (14) and (15) and using Eqs. (13), we obtain 



dn 



12 



902 



+ 



2V d(t>i ^ 3(y2 + ^i) 



901 



Analogously, we derive the expression for dn22/d^: 

dn22 ^ 1 <902 _ 2QV 90i ?,{QV^ + 62) 

- 302 (Q\/2 _ 3^2)2 dr (QV - 362)^ ' 

Differentiating Eq. (18) over ^ and using expressions (19) and (20), we obtain the KdV equation. 



(14) 

(15) 
(16) 

(17) 
(18) 

(19) 



(20) 



dr 



- + ^01 ^r— + -^zz:r = 0, 
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(21) 



where 



3fii{V^ + 9i) 3fi2{QV^ + 02) 
(y2 _ 3^^)3 (gy2 _ 3^2)3 



- 1 



-7 



/^i 



+ 



QAi2 



After switching to the coordinate r] = ^ — ut, the solution to Eq. (21) that satisfies the condition 
01 (^ ±00) = can be rewritten as 



01 = 0^ch-'((M/26)^/2^), 0^ = 3m/A 



(22) 



For numerical analysis, it is convenient to introduce the dimensionless coefficients a — /i2//^i = 
ni2o/niio and P — Znao/nao, which are the density ratio between the negative and positive ions 
and the ratio of the dust grain charge to the charge of a positive ion, respectively. Figure 1 shows 
how the wave amplitude depends on a at different values of /3. It is seen that the critical density 
of negative ions decreases in the presence of a negatively charged dust. In contrast, if the dust is 
charged positively, then the critical density increases. If the dust is charged negatively, then, for a 
given value of the parameter a, the soliton amplitude in a dusty plasma is higher than that in the 
absence of dust, and vice versa, if the dust is charged positively, then the soliton amplitude in a 
dusty plasma is lower than that in a bi-ion plasma without dust. 



4 MODIFIED KORTEVEG-DE-VRIES EQUATION 



Let us now assume that the density ratio between the negative and positive ions is equal to its 
critical value, i.e., A — 0. In this case, in order to describe the nonlinear properties of the medium, 
it is necessary to take into account the higher orders of the expansion in the small parameter. We 
again apply the method of coordinate stretching, 

^^e{x-Vt), T^eH. (23) 

In the first order in £, the equations remain unchanged. In the second order in £, the d/dr derivative 
in Eqs. (14)-(17) is absent. This allows us to integrate relationships (19) and (20) analytically. Using 
continuity equations (14) and (16), we also obtain expressions for V21 n V22 that are necessary to 
construct the theory in the third order in e: 



^2-3^1 2(\/2_3^^)3<^l' ^2_3^/2,+2(^2_3^^)3<^l' 



02 , S{QV' + 92) 2 V V{QV' + 992) 2 (^,. 

^22 Qy2 _ 3^^ + 2(Q\/2 _ 3^2)3*^" QV^ - 3^2 2(gV2 _ 3^2)3'^'' ^ ^ 

In the third order in e, Eqs. (l)-(4) take the form 

drill ^.driis dv^ d 

-g;^-V-^ + -^ + —{ni2Vii+niiVi2) = 0, (25) 

dvn avi3 d N , "^^s ni3 <9 

- + ^(^11^12) + ^ + + 3^i^Ki^i2) = 0, (26) 

dn2i ^.dn23 . dv23 ^ d 

-g;^-V-^ + -^ + —{n22V2i+n2iV22) = 0, (27) 

9^21 ..9^23 . d 1 d(p3 3^2^23 3^2 «9. . . . 

Excluding the terms with dvi^/d^ h dv23/d^, we obtain 

dnis ^ 1 d(l>s 2V d(t>i 3(y2 + g,) d 

V - Wi {V - 3^1)2 dT (y2 _ 3^,)3 ^'^i'^2j+ 

+ 2(1/2 _ 3^^)5 'Pi^- (29) 
dn^ ^ 1 >903 _ 2gK g0i 3(gi^2^g^) g 

9^ gi/2 _ 3^2 9e (gv^2 _ 3^2)2 dr (g\/2 _ 3^2)3 ^"p^^^) 

3{5Q^V^ + 30QV^92 + 9el) ^,d<Pi 

2{QV' - 362)' dC ■ ^ ^ 

Differentiating Poisson's equation over $, in the third order in e yields 



^-p- = 03 + 0102 + 0i/ 6 - Hini3 + /i2?^23- 



(31) 



Then, substituting expressions (29) and (30) into Eq. (31), we arrive at MKdV equation 



(32) 



where 



3/ii(51/^ + 30^2^1 + 961) 3/X2(5QV^ + 30QV^92 + 99^) 



+ 



Switching to the coordinate r] = ^—ut and integrating Eq. (32) under the condition 0i(^ ±cxd) = 0, 
we obtain a solution in the form of a solitary wave. 



01 = ±0„ch-i((2M/6)i/2^), 0^ = {Qu/B) 



1/2 



(33) 



Let us now assume that the negative ion density is close to its critical value but is not exactly 
equal to it. Then, the second-order terms in Poisson's equation yield 



ne2 - piinu + A*2n22 = 02 + 0i/2 - Ail 



+ 



+A*2 



+ 



-A<Pl. 



(34) 



_QV^-3e2 ' 2{QV^-3e2f^\ 
If A = 0(e), then we have A(j)1 — 0{e^). Accordingly, there is no terms of the second-order in e on 
both the left-hand and right-hand sides of Poisson's equation, because, although each of the terms 
02, 01, ni2 and ^22 is on the order of e^, their sum is a quantity of the third order in e. In this case, 
instead of Eq. (32), we obtain the equation 



901 , 901 D,2 901 , ^•9^01 



+ 77- 



0, 



(35) 



dr 9^ 2 

where the second and the third terms on the left-hand side are of the same order (note that, in 
deriving this equation, the term A ■ 9(0i02)/9^ was omitted because it is of the higher order than 
the other terms). 

In the terms of the coordinate r] = ^ — ur , the solution to Eq. (35) that has the form of a solitary 
wave and satisfies the condition 0i(^ ±00) = is expressed as 



6ti/A 



1 ± Vl + 6Bu/A'^ch[{2u/b)y^r]] ' 



The amplitude of the soliton is 



6u/A 



(36) 



(37) 



l±y/l + 6Bu/A^' 

If S = or A = , then expressions (36) and (37) transform into formula (22) or (33), respectively. 

As was mentioned above, the critical value of the negative ions density is unambiguously related 
to the dust density. The value of /3 at which the nonlinear coefficient A in Eq. (21) is zero is plotted 



in Fig. 2 as a function of the parameter a for different mass ratios between negative and positive 
ions. The points at which the vertical line /3 intersects the plots correspond to the critical densities 
of negative ions in a plasma without dust. Figure 3 shows the amplitude of an MKdV soliton as a 
function of /? at different values of Q. A numerical analysis shows that the coefficient B is positive 
at any value of Q. Therefore, in the vicinity oi A — 0, Eq. (35) describes solitary waves and has no 
solutions in the form of a double layer. 



5 ANALYSIS OF THE QUASI-POTENTIAL 



Let us consider the propagation of a wave of arbitrary amplitude. Passing to the frame of reference 
related to the wave, r) — x — Mt (where M is the Mach number in terms of c^); integrating Eqs. 
(l)-(4) and substituting the resulting expressions for the ions densities into Poisson's equation, we 
obtain 

2 



C/(0) = 1 - e-^ + 



(38) 



1+ 



^2f^i+ 



a, 



3 

2- 



'2+ 



1 + ^) 



3/2 



6V392 

where (Ti± = M ± ^/3^, a2± = VQM ± y/W^. 

Expanding the expression ^7(0) in a power series in yields 



U(6) ^ --( I- ^2 

2 V M2 - 3^1 QM^ - 3^2 



+ 



0V 



+- 



1 /3/ii(M2 + Oi) 3/i2(gM2 + 62) 



6 V (M2 - 3^i)3 (gM2 - 3^2)^ 



(39) 



1 

+ 24 



3l^i{5M'^ + 30i\Pei + 991) 3/X2(5g2M^ + 30QM^e 



+ 



1 



(M2 - 3^1)5 (gM2 - 3^2)^ 

In deriving nonlinear equations (21) and (32) (see Sections 3 and 4), we set ^ = a; — Vt, where 
V is the propagation velocity of a linear wave. To find weakly nonlinear solutions, we then passed to 
the variable rj — x — ut, where is a nonlinear correction to the propagation velocity. Therefore, for 
the subsequent analysis to be consistent with our previous results, it is necessary to set M = y + u, 
where u ^ V . In this case, the expression in the parenthesis in front of 0^ equals u/b and, in the 
other terms, it is sufficient to replace M with V. As a result, we obtain: 



driJ 3b^ 3b^ 



(40) 
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The solution to this equation coincides with expression (36). 

If the wave properties are determined by the coefficient A and the nonlinear correction u to the 
propagation velocity is not too small, then the higher order in the expansion of the Mach number 
in the power series in u/V should be retained in the coefficients of expression (39). In this case, the 
solution to Eq. (38) can be written as 

= f ch-2(r;/A) + ^{^EbV - B/Ay ch-'{7^/A) - ch-2(r;/A) th'{7^/A), (41) 

where A = {2b/uy/^{l + SAbu/A) and 

p + 3^1 ^ Qv^ + 3e2 . V + e, ^ n + ^2 

-^ = /^i7T7^ TTTT^I - ^^Q 77^ ^TT^' ^ = 17775 ^TTT^^ + ^^Q: 



Note that taking into account the next term in the expansion of M is equivalent to constructing 
a KdV equation that is nonsecular in the second order in e and has solutions in the form of the 
so-called „dressed" solitons [13, 14]. From the standpoint of such solutions, the first term in Eq. (41) 
correspond to the „core", whereas the second and the third ones form the „ cloud". 

For numerical analysis, it is convenient to use the quantity M' — M/V (the Mach number divided 
by linear propagation velocity (13), which is a root of a biquadratic equation). Far from the critical 
density of the total negative charge of the dust and ions, there is only one type of solitary wave — either 
compression or rarefaction one. At the critical density, both types of wave can exist simultaneously 
and the plot of the quasi-potential U{(f)) intersects the axis @[phi] at the points that are symmetric 
relative to = 0. When M' is close to unity, the amplitudes of these waves are small and are well 
described by formula (33). As M' ncreases, the amplitudes of both positive and negative waves of 
the potential increase. Figure 4 illustrates how the shape of the quasi-potential well for the critical 
values a = 0.202 and P = 0.2 (all the other parameters are specified in the figure caption) varies 
with increasing Mach number. As the Mach number increases, the wave amplitude increases and, at a 
certain value M^^ (for the parameters corresponding to Fig. 4, it is M^^ Ki 1.129) , the kinetic energy 
of an incoming ion in the soliton frame of reference becomes less than its potential energy at the 
center of the soliton; i.e., the ion is reflected. Figure 5 shows the shape of the quasi-potential well for 
the case where the dust density (or the negative ion density) is not equal to its critical value. As an 
example, the figure illustrates the situation in which the negative ion density is constant, whereas the 
dust density increases and the wave transforms into a rarefaction soliton {(p^ < 0). Slight variations 
in /3 leads to insignificant variations in the amplitudes of compression and rarefaction solitons. For 
small-amplitude solitons, this is clearly seen from formula (39). If ^ 6Bu, then the soliton 
amplitude is (f)m ~ ±(6m/5)^/^(1 =f A/ \J 6Bu). In the case at hand, we have ^4 < 0; therefore, the 
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amplitude of the compression soliton is slightly higher than that of the rarefaction one. It can be 
seen from the figure that, for waves of arbitrary amplitude, the difference between the amplitudes 
of the compression and rarefaction solitons becomes considerable. As the dust density increases, the 
amplitude of a positive-potential soliton becomes higher than i.e., the condition for the 

existence of such solitons is violated. In this case, the quasi-potential has only one turning point, 
which corresponds to a rarefaction wave. As the dust grain charge decreases below its critical value, 
the amplitude of a negative-potential soliton increases and the turning point of the quasi-potential 
for such solitons also disappears when /? decreases below a certain value. Xote that the range of the 
values of the dust-density parameter (3 within which both types of solitary waves exist simultaneously 
is rather wide. Thus, in the case at hand, the parameter (5 can vary in the range 0.159 — 0.42. If (3 
is lower than —0.159 (which corresponds to a positively charged dust), then only positive-potential 
solitons can propagate in the plasma. In contrast, if /3 is higher than 0.42, there are only solutions 
in the form of rarefaction solitons. The table presents the ranges of f5 within which compression and 
rarefaction solitons can exist simultaneously at different values of a. 

It should be noted that all the numerical calculations of the wave amplitudes and nonlinear 
coefficients were performed for the fast ion-acoustic mode [6]. An analysis of the numeric results has 
shown that the phase velocity of slow ion-acoustic waves is on the order of the ion thermal velocity; 
therefore, these waves undergo strong Landau damping. 

6 CONCLUSIONS 

Thus, using the approach based on the KdV equation, we have shown that the presence of a 
negatively charged dust in plasma leads to a decrease in the critical density of negative ions. For a 
given value of a (the ratio between negative and positive ion densities), the amplitude of a compression 
soliton A > 0) in a dusty plasma is larger than that in a bi-ion plasma without dust, whereas the 
amplitude of a rarefaction soliton [A < 0) is smaller. In contrast, if the dust is charged positively, 
then the amplitude of compression (rarefaction) solitons in a dusty plasma is smaller (larger) than 
that in a plasma without dust. 

When the dust is charged positively, there always exists a critical value of the parameter a (the 
density ratio between the negative and positive ions). For a negatively charged dust, the condition 
A — can be satisfied for /3 values below a certain critical value Per, at which a becomes zero and 
above which only rarefaction solitons can propagate in the plasma. In the absence of negative ions, 
the critical value Per in a cold plasma is equal to 2/3 [10] and it increases slightly if the thermal 
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effects are taken into account. MKdV solitons in a dusty plasma exist along the corresponding curves 
in Fig. 2 and only for (5 < (5cr- 

An analysis of the quasi-potential shows that, for waves of arbitrary amplitude, the curve that 
separates the domains of existence of compression and rarefaction solitons in the plane /9 — a is rather 
conventional and that there is a broad transition band on both sides of this curve within which both 
positive- and negative-potential solitons can exist simultaneously. Below this curve (in the domain 
of existence of compression KdV solitons), the amplitude of rarefaction solitons is larger than that 
of compression ones. On the other hand, above this curve (in the domain of existence of rarefaction 
solitons), the amplitude of compression solitons is larger. Above the transition band, only rarefaction 
solitons exist, whereas below this band, only compression solitons can propagate in the plasma. 
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ymax 
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n 1 
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n ot? 
O.oob 




0.3 


-0.902 


0.251 




0.3 


-0.077 


0.400 




0.3 


-0.240 


0.436 




0.4 


-2.204 


0.078 




0.4 


-0.290 


0.290 




0.4 


-0.408 


0.335 




0.5 


-4.113 


-0.115 




0.5 


-0.520 


0.177 




0.5 


-0.578 


0.232 



Ranges of the parameter /3 within which both compression and rarefaction solitons can exist for 
Oi = 02 = 0.01, M' = 1.05, and different values of Q and a. 
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FIGURE CAPTIONS 



Fig. 1. Amplitudes (prn of (a) compression {A > 0) and (b) rarefaction solitons {A < 0) as functions 
of the parameter a — ni2o/?T'iio for Q — 32, 9i — 02 — 0.05, u — 0.01 , and different values of /3. 

Fig. 2. Curves corresponding to the condition ^4 = in the plane {/3 — a) for 9i — 02 — 0.05, 
u = 0.01, and different mass ratios between the negative and positive ions. 

Fig. 3. Amplitude of an MKdV soliton vs. (3 for 6i = 62 = 0.05, u = 0.01, and different values of 

Q. 

Fig. 4. Shape of the quasi-potential well at the critical density of negative ions for a — 0.202, 
/3 = 0.2, Q ^1,01 = 92^ 0.01 , and different values of the Mach number M: 1-1.03,2-1.05,3-1.1. 

Fig. 5. Shape of the quasi-potential well for a = 0.2, Q = 1, 9i = 62 = 0.01, M = 1.05, and 
different values oi (3: 1 - 0.15, 2 - 0.25, 3 - 0.35. 
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Fig 2. 
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